Fractional Reaction-Diffusion Equation 



Kazuhiko SekiQ Mariusz WojcikQ and M. Tachiy£0 

National Institute of Advanced Industrial Science and Technology (AIST) 
AIST Tsukuba Central 5, Higashi 1-1-1, Tsukuba, Ibaraki, Japan, 305-8565 
(Dated: February 2, 2008) 

A fractional reaction-diffusion equation is derived from a continuous time random walk model 
s ! ' when the transport is dispersive. The exit from the encounter distance, which is described by the 

algebraic waiting time distribution of jump motion, interferes with the reaction at the encounter 
distance. Therefore, the reaction term has a memory effect. The derived equation is applied to the 
geminate recombination problem. The recombination is shown to depend on the intrinsic reaction 
rate, in contrast with the results of Sung et al. [J. Chem. Phys. 116, 2338 (2002)], which were 
. obtained from the fractional reaction-diffusion equation where the diffusion term has a memory 

' effect but the reaction term does not. The reactivity dependence of the recombination probability 

is confirmed by numerical simulations. 
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I. INTRODUCTION 



Anomalous diffusion processes occur in many physical systems for various reasons including disorder in terms of 
energy or space or both 1,2 . The anomalous diffusion processes are expressed as (r 2 (t)) oc t a for the displacement 
r(t); the process is called sub-diffusive or dispersive if < a < 1 while it is called super-diffusive if a > l-' 2 . One 
d . of the useful theories for such processes is the continuous time random walk where the non-Poissonian waiting time 
distribution of jump motion is introduced^. For sub-diffusion processes the fractional diffusion equation has been 
derived from the continuous time random walk model with the power-law waiting time distribution function™. The 
theory is continuous in space and thus is very useful to introduce the effect of interactions and the boundary conditions. 

The fractional diffusion equation is valuable for describing reactions in the dispersive transport medial. In 
the theory of diffusion-controlled reactions the boundary conditions describing reactions are very important. For 
the perfectly absorbing boundary conditions the generalization of the Smoluchowski theory based on the ordinary 
diffusion equation to the fractional diffusion equation is straightforward. On the other hand, for the partially 
reflecting boundary conditions the exit from the encounter distance interferes with the reaction at the encounter 
distance, therefore the latter is influenced by the non-exponential waiting time distribution function. Thus, the 
generalization of the conventional partially reflecting boundary conditions associated with the ordinary diffusion 
\Q , equation to those associated with fractional diffusion equation is not obvious. Recently, Sung et al. phenomeno- 
logically introduced the fractional reaction-diffusion equation and found that the recombination probability of a 
particle starting from ro for the partially reflecting boundary condition is equal to that obtained for the perfectly 
absorbing boundary condition 6 , and thus independent of the intrinsic reaction rate at the encounter distance. When 
the partially reflecting boundary condition is used, a fraction of particles that arrive at the encounter distance will 
recombine, while the others will escape the reaction. Therefore, the recombination probability for the partially 
reflecting boundary conditions should be lower than that for the perfectly absorbing boundary condition. Fractional 
reaction-diffusion equations or continuous time random walk models are also introduced for the description of 
nonlinear reactions, propagating fronts and two species reactions in sub-diffusive transport media 7 . However, 
these fractional reaction-diffusion equations share the same structure as that used by Sung et al. ; the reaction 
term does not have any delay effect and only the diffusion term has a memory effect. In this paper we derive a 
fractional reaction-diffusion equation from a continuous time random walk model ; it gives the proper partially 
reflecting boundary conditions applicable for the fractional diffusion equation. The reaction term as well as the 
diffusion term has a memory effect because the exit from the encounter distance, which is described by the algebraic 
waiting time distribution of jump motion interferes with the reaction at the encounter distance. Our theory is 
different from the more macroscopic description of nonlinear reaction in sub-diffusive transport media™, where 
the reaction is accounted for not by a space dependent sink term but by the law of mass action which is space 
independent. The recombination probability is obtained from the fractional reaction-diffusion equation thus de- 
rived. It is also shown by simulations that the recombination probability indeed depends on the intrinsic reaction rate. 
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II. FRACTIONAL REACTION-DIFFUSION EQUATION 



We consider geminate recombination of a B particle starting at f*o with A. B particle migrates toward A by 
anomalously slow diffusion, (r 2 (t)) oc t a with a < 1. Reaction takes place at the encounter distance of R to R + dr. 
After B particle leaves the encounter distance it performs random walk described by the sub-diffusion kinetics. When 
the migration dynamics is described by the ordinary diffusion equation, the reaction is accounted for by assuming that 
the flux of B into the region R to R + dr is proportional to the density of B at R and the resultant boundary condition 
is well established^. This boundary condition is referred to as the partially reflecting or radiation boundary condition. 
It is also well known that imposing the boundary condition is equivalent to introducing the sink term in the diffusion 
equation (reaction-diffusion equation) with the perfectly reflecting boundary condition^. However, if the motion of B 
particle is sub-diffusive, the corresponding boundary condition is not known, though the sub-diffusion motion itself 
is described by the fractional diffusion equation. On the other hand, in the lattice model the sub-diffusion motion is 
derived from the theory of continuous time random walk and reaction can be easily accounted for in this theory. We 
start from the lattice model and take the continuous limit in order to obtain the sink term for the fractional diffusion 
equation. 

The continuous time random walk is specified by the waiting time distribution, ip(t), of making a jump to a 
neighboring site at time t in the absence of reaction. To be more specific a particular model is introduced where the 
algebraic asymptotic tails result from the distribution of the site energy^. The power-law waiting time distribution 
is also obtained from an exponential distribution of inter-trap distances together with an exponential dependence of 
the jump rate on the inter-trap distance^. As long as the long time behavior is concerned our conclusion is model- 
independent. One of the site energy distribution functions which we encounter most frequently is the exponential 
distribution 9 , 

g(E) = cxp [-E/(k B T c )\ /(AfeTo). (2.1) 

For the activated release rate, 

1 (E)= lt exp[-E/(k B T)}, (2.2) 
the waiting time distribution for release is given by 9 - 

m = r dEg{E) 1 {E)e^{- 1 {E)t) ~ ^f" + * } , (2.3) 
Jo It 1 

where a = T/T c and T(z) is the Gamma function. Since we are interested in dispersive transport a < 1 is assumed. 
In the small s limit the Laplace transform of the waiting time distribution function can be expressed as^ ^(s) ~ 
1 — [ira/ sin(7ra)] (s/j r ) for s/j r < 1. Two types of waiting time distribution are defined at the encounter distance 
; one is the waiting time distribution function of making a jump to a neighboring lattice site ip ou t(t) and the other 
is the waiting time distribution function of reaction ij) TC {t). The waiting time distribution of making a jump at the 
encounter distance is given by, 

^out(i) =^(t)exp(-7rct). (2.4) 

The waiting time distribution of reaction is defined as the probability that the particle which is initially at a site 
in the encounter distance will undergo reaction without making a jump at time t. It is given by the reaction rate 
constant, 7 rc , multiplied by the remaining probability of particles at the site in the reaction zone, which decays either 
by jump motion or reaction, 

dtxip(h). (2.5) 

It is implicitly assumed in our model that in the presence of site energy distribution g{E) the reaction takes place 
from any energy level with the same rate 7 rc . 

We define the vector characterizing a jump to the nearest neighbor site j, by bj (j = 1, 2, • • ■ , 2d) and the jump 
length b. In the limit of small A1Z the region from R — A1Z/2 to R + A1Z/2 can be regarded as the encounter distance. 
In this section we first consider the system without the reflecting boundary condition and introduce later the reflective 
sphere of radius r = R. The equation governing the probability rj(fi,t) of just arriving at site r*j at time t can be 
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written as, 
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2d- 



/ dh ^(t-tiWn-bjih)- 
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<L r - Slp(t 



h)AH(R - AK/2 - r k ; AK^f^h) 



+ 6? tA> 6(t), (2.6) 

where AH(R — r; A1Z) = [H(R + A1Z — r) — H(R — r)] and H{r) is the Heaviside step function, namely, H(r) = 1 
for r > 0, otherwise H(r) — 0. Sip(t) = ip(t) — ip ou t(t) is the difference between the waiting time distribution of 
making a jump at other sites and that at the encounter distance. The difference arises because of the reaction at the 
encounter distance. The probability of just leaving site fi at time t is given by f Q dt\i]){t — ti)ij(fi, t\). By subtracting 
this quantity from both sides of Eq. I|2.6|) we obtain the balance equation. After introducing the Laplace transform, 
■0(s) = Jq 00 dtexp(—st)ijj(t) 7 the balance equation is written as, 



2il 



2^£^Wi ~bj,s) -fj(fi) 



^ 2d 

^{s)-: J2 AH ( R - A7Z / 2 - r> - h ; An) f)(fi - b j)S ) + S r ^ . 



2d 

i=i 

Now we introduce the perfectly reflective sphere of radius R. Reaction takes place when a particle enters the reaction 
zone, which is defined by the volume between R and R + A1Z/2. The shortest distance between a point outside 
the reaction zone and the reflective sphere is A1Z/2. Therefore any particle which makes a jump from outside the 
reaction zone and is reflected by the sphere cannot reach a point outside the reaction zone, unless the jump length is 
larger than All. Thus, we set A7L = b to guarantee that any particle which makes a jump from outside through the 
boundary of the reaction zone lands within the reaction zone. As we will show later by simulations, this choice of the 
reaction zone also satisfies the condition that practically all the particles which make a jump from inside the reaction 
zone land outside it. Therefore, the particle in the reaction zone comes from outside it, 



1 - ip(s) fi(ri,s) = ip(s) 



^ 2d 

2d^ f]{ 



n - bj, s) - fi(fi) 



Sft f for R < ri < R + - , 



(2.7) 



where fi — bj stands for the position which leads to fi after making a jump bj. If \fi — bj\ > R, the particle is not 
reflected by the sphere of radius R, therefore fi — bj equals to fi — bj when \fi — bj\ > R. On the other hand, if 
\fi — bj\ < i?, the particle is reflected by the sphere. In this case fi — bj is given by the position shown in Fig. 1 a) 
and b). Strictly speaking, even if \fi —bj\ > R, there is a case where a particle is actually reflected by the sphere, as 

bj. But such events are obviously rare. The balance equation outside 



shown in Fig. 1 c). In this case fi 
the reaction zone obeys, 



bj ^ n 



1 - ij>(s) fj(fi,s) = 



1 
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3=1 



bj,s) -fj(fi) 



2d 
3=1 



R 



; 2 ) v(n 



b 3 ,s) + 5? i: p 



for R + - < n 



(2.8) 



We introduce the corresponding probability density, r\ (f, t) — J^. rj {fi,t) S (r ' — fi) . Eq. I|2.7f) does not directly involve 
the effect of reaction, while Eq. (|2.8I) does. Moreover, for a small jump length the contribution from Eq. (|2.7|) is 
negligible. Therefore, we focus our attention on Eq. Ij2.8|l . Performing Taylor expansion, the leading terms satisfy 
the differential equation, 



1 — ip(s) r)(r, s) 
b 



-■0(s)V 2 7} (r, s) 



5?p(s)5 (R — r)f] (f, s) + S (f— ro) , 



(2.9) 
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where the following approximation is employed, 

1 2d f 



3=1 
1 



f — 6, 



1 2d 



f — b. 



3=1 



2 AH(R + b/2 - r; b) = ^6{R + b/2 - r) ~ ^S(R - r). 



(2.10) 



The first equality follows from the definition of r — bj, which is equal to r — bj when R < r — bj < R + b/2, while it 



is given by the corresponding reflected vector when R — b/2 < 



< R as described in Fig. 1 b). The quantity 



in the summation becomes AH (R + b/2 — r; b) for bj directed outward from the sphere of radius R and vanishes 

otherwise, and among 2d of bj vectors half of them are directed outward. This explains the second equality of Eq. 
(|2.1U|) . The third equality is due to an approximate expression of the delta function given by the first derivative 
of the Heaviside step function. So far, we have considered the probability density of just arriving at f at time t 




R-b/2 R R+b/2 



FIG. 1: The position r — bj which leads to f% after making a jump bj. a) for |f — bj\ < R and R < r < R + b/2. b) for 
\f — bj\ < R and R + b/2 < r. c) for \r — bj\ > R. In this case, although \r — bj\ > R, the particle is reflected and f — bj is 
actually not equal to r — bj . 

, r)(r,t). The usual probability density is defined in terms of the remaining probability, cf>(t) = J. dtxip(tx) as 
p (F, t) = J * dt\4> (t — ti) rj (r, ti) . By noticing the relation, cj>(s) = (l — ?/>(s)^ / s, Eq. (|2.9() becomes 

sp(r, s) - p (f, t = 0) 



l-^(s) 

^( S )V 2 /5(f, s) - 2TrR 2 b6i[,( S ) S{7 i ~ R f ) P {r, s)\ . (2.11) 

Since a particle cannot penetrate the sphere of radius R, the perfectly reflecting boundary is imposed at r = R. In 
the long time limit which is expressed in the Laplace domain as s ~ + we have 1 — %/>(s) — > s a and Eq. (|2.11ll reduces 
to 

sp(r, s) - p(r f ,t = 0) 



&2 V 2 p(f, a) - 2nR 2 b^ c (0) 5{r A „f } p(r, s) ) . (2.12) 



l-V(a) 

2 ^,,. ^^ ! ^c(0)^ 2 

where V'rc(O) = 1 — - 0out(O) is introduced, which follows from the fact that particles at a given site in the reaction zone 
perform either jump or reaction, namely ip IC (t)dt + ipout{t)dt = 1 ; the overall probability of reaction, ■0 rc (O), 
is the branching ratio of undergoing reaction at a site in the reaction zone and V'out(O) represents the branching ratio 
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of making a jump at the same site. For the normal diffusion, a = 1, Eq. I|2.12[l yields the usual reaction-diffusion 
equation after inverse Laplace transform. For sub-diffusive transport media, a < 1, the right-hand side of Eq. (|2.12|) 
is given in the time domain by time differentiation of the convolution of the retardation function of algebraic form 
with the function which includes both the diffusion term and the reaction term. Now, we investigate this feature. 

The Laplace transform of V'rc(i) defined in Eq. (|2.5|l is obtained as ^> rc (s) ~ 1 — ip (s + 7 rc ) . In the limit of fast 
release rate 7 r and fast reaction rate 7 rc we get 

i c ( S )~^-(^y (2.i3) 



sin 7r<y \ 7, 

for reaction-limited condition, 7 rc < j r . In this case ip T c(s) is a constant independent of s. On the basis of the 
diffusion-reaction model described above we define the generalized diffusion constant D a as 



Then the inverse Laplace transform of Eq. I|2.12|l is obtained as 

9 , s 9 f* , 1 

di pir > t) = aij dtl nc 



D a V 2 p (r, tx) - k « S ^ nR ? P (r, t x ) 

where the generalized intrinsic reaction rate k a is defined as , 

k a = ~f? c 2nR 2 b. (2.16) 

The exit from the encounter distance, which is described by the algebraic waiting time distribution of jump motion 
interferes with the reaction at the encounter distance. Therefore, the reaction term has a memory effect like the term 
describing sub-diffusion motion. So far b has been assumed to be small but finite and we have investigated equations 
under the condition of large 7 r and 7 rc . In order to derive the equations valid in the continuous limit we have to take 
the limit of b — > 0, j r — > oo and 7 rc —> oo with D a and k a kept finite. Here, the fractional reaction-diffusion equation 
has been derived for a specific initial condition, but the linearity of the model guarantees a much wider range of 
applicability: essentially, any initial distribution will be acceptable. 

In the reaction-diffusion equation, Eq. I|2.15|l . the reaction is accounted for by introducing the sink term, with the 
perfectly reflecting boundary condition imposed. If the reaction is accounted for by imposing the boundary condition 
associated with the fractional diffusion equation!, 

m" (r ' t} = I L * dtl ffe 7rh^ DaV2p (r ' h) ' (2 - 17) 



){t-t 1 f- a 

(2.15) 



what is the relevant boundary condition? By multiplying both sides of Eq. (|2.15|) by J fl e Airr 2 dr we obtain in the 
limit of a small e, 



d 

D a -^p{r,t) 



% P {R,t). (2.18) 



Therefore the boundary condition for the fractional diffusion equation, Eq. 1)2. 17)) . is the e — > limit of Eq. I|2.18|l . 
Eq. (|2.15() or Eq. (|2.17|l with Eq. (|2.18() is the most important result of this paper. 



III. RECOMBINATION PROBABILITY 



Now we apply our equation to the geminate recombination problem. The recombination probability is obtained by 
the usual procedure. The fractional diffusion equation can be cast into the form, -^p(r,t) = —divj(r,t), where the 

rate of the change of the density is related to the current defined by, j(r, t) = — Jj J ^^lj^y ( t _ t] y- a D a gia,Ap (r, t\) . 
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The recombination probability k(?*o) is expressed as At(ro) = - J °° dtiA-KR 2 j{R,ti) ■ R/R, which can be calculated 
by the standard method as, 

The familiar form of the recombination probability is derived, which clearly shows its dependence on the intrinsic 
reaction rate, in contrast with the results of Sung et al£, which were obtained from the fractional reaction-diffusion 
equation where only the diffusion term has a memory effect and the reaction term does not. One should also note 
that both the diffusion constant and the intrinsic reaction rate scale with parameter a as shown in Eqs. (|2.14() and 
(|2.1tj|) . The scaling for the intrinsic reaction rate is due to the reaction model used. For other models different scaling 
with a may be derived. 

Eq. p.lfl can be derived more directly in the following way. Since the probability that a particle which starts at ro 
will visit the spherical shell of radius R is equal to the recombination probability for the perfectly absorbing boundary 
condition^, it is given by R/rQ. Inside the reaction zone the probability that a particle makes a jump without reaction 
is V'out(O). After a jump the particles in the reaction zone may still remain inside the reaction zone or leave the 
reaction zone. In the limit that the width of the reaction zone goes to zero, we consider two probabilities, namely, the 
probability that a particle at the spherical shell of radius R will make a jump to that of radius R + b with b being the 
jump length, mip ou t(0), and that of making a jump to another position at the reaction distance R, mp ou t{0), where 
m and n are the branching ratios of respective jumps and m + n = 1. On the other hand, the probability that it will 
undergo reaction is "0 rc (O). Some particles may recombine at the first visit to the encounter distance R. Other part 
of particles escape recombination at the first encounter and make a jump to the sphere of radius R + b or make a 
jump to another position at the encounter distance. In the former case some particles may recombine at the second 
encounter after jumping from the sphere of radius R + b to the encounter distance or again escape recombination at 
the second encounter. In the latter case particles may recombine at the same position or make another jump, and so 
forth. Accordingly, the probability that a particle which starts at ro will ultimately undergo reaction is, 

R f R * \ R 

t(r ) = ^rc(O)— + ^rc(O) — -rTOVw(O) + m/w(0) )—+••■ 

ro \R + b J r 

R ^ rc (0) 



r °l- (mtjj out (Q) R + Ww(0) 



R 1 



r R + nb mb 



R + b R + bijj ic (0) 
R 1 



ro 



(3.2) 



1 



R AM 

Therefore, the recombination probability is given by Eq. I|3.1|l with k a generalized to, 

k a ~ 2TrR 2 b 1 */m. (3.3) 

In simulations we cannot take the zero limit for the width of the reaction zone and the value of m depends on the 
width. By comparing to the simulation results we find m = 1 for the width of 6/2 and m = 1/2 for the width of 6, 
as will be shown later when we compare the analytical result with simulations. Eq. Ij2.15|) is derived for the width 
of the reaction zone equal to 6/2 before taking the limit that the jump length goes to zero. Since m = 1 is obtained 
for this choice of the width, the intrinsic reaction rate Eq. (|2.16|) is consistent with the definition of the generalized 
intrinsic reaction rate, Eq. (|3.3(l . 



IV. SIMULATIONS 



In order to confirm the validity of the above results we have carried out simulations. In our calculations the random 
walk is realized as a sequence of instantaneous jumps between the traps, with the trap energies generated according 
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to the distribution Eq. i|2.1|) . The rate of release from the traps is assumed in the form given by Eq. (|2.2|l , and the 
detrapping time for a trap with energy E is obtained from the exponential distribution with the mean value l/"/(E) . 
We assume the Gaussian distribution of jump lengths and calculate displacements of particle B along each Cartesian 
coordinate as Ax, Ay, Az = b x A^(0, 1), where N(Q, 1) represents a random number obtained from the standard 
normal distribution. The trajectory of B is calculated until it either reacts with A or escapes to a large distance r max 
. The simulation is repeated for a large number of independent trajectories, which allows us to calculate the reaction 
(and escape) probability. 

We assume that the reaction occurs within a thin spherical shell R — AR/2 < r < R + AR/2 , where R is the 
reaction radius, and is characterized by the reaction rate 7 rc . In the simulation the reaction is modeled in the following 
way. When a particle B is trapped within the reaction shell, we generate the reaction time t IC from the exponential 
distribution with the mean value l/"f rc , and compare it with the detrapping time t calculated for the current trap. 
If < rc < t then the reaction occurs, otherwise the particle jumps to another trap and the simulation is continued. 
A reflective boundary is established at r = R — AR/2 , and the particle is bounced back when a jump across this 
boundary is attempted. The reaction model described above represents a partially diffusion-controlled reaction with 
the classical second-order rate constant k — 47ri? 2 Ai?7 rc when the transport is described by the normal diffusion. 

The calculations presented in this work were carried out for ro — 2R and r max = 100i?. Error in the calculated 
escape probability due to the finite value of r max is estimated as about 1 %. The number of independent trajectories 
generated in each simulation was at least 10 . The jump length parameter was assumed as b = 0.17?, and the values 
AR — b and AR = b/2 were used. In the calculations we used reduced units, with R taken as the unit of length and 
l/7r as the unit of time. The calculated diffusion constant is in good agreement with Eq. I|2.14[l . The simulation 
results also show that the escape probability 1 — k (r ) depends on the intrinsic reactivity not only for the normal 
diffusion, but also in the sub-diffusive case as shown in Fig. 2. The analytical result of Eq. i|3.1fl coincides very well 




FIG. 2: The escape probability versus log 10 {RD a /k a ). Initial distance is tq = 2R. The symbols denote simulation data. 
The closed circles, closed triangles and pluses correspond to the a values of 0.5, 0.75 and 1.0, respectively for the width of the 
reaction zone b/2. The open circles, open triangles and crosses correspond to the a values of 0.5, 0.75 and 1.0, respectively for 
the width of the reaction zone 6. The line shows the analytic result of Eq. 13.11 . 



with the simulation data. When the width of the reaction zone is b/2 we use m = 1 and when the width of the 
reaction zone is b we use to = 1/2 for the definition of the generalized intrinsic reaction rate, Eq. (|3.3|) . Since the 
simulation data coincide very well with the analytical result of Eq. I|3.1|l the above choice of m is also justified. When 
the width is the same as the jump length of the particle motion, b, the particles which make a jump from inside the 
reaction zone and are reflected by the sphere of radius R still remain inside the reaction zone. Approximately half 
of the particles make a jump outward from the sphere of radius R and leave the reaction zone after a jump. On the 
other hand, when the width is as small as b/2, almost all the particles including those reflected by the sphere of 
radius R leave the reaction zone after a jump. 
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V. RECOMBINATION PROBABILITY IN A LATTICE MODEL 



Finally we analyze the reaction kinetics in the full lattice model where the reaction proceeds only at the origin. 
The recombination probability of a particle starting from r*o for the partially reflecting boundary condition is obtained 
from the known result for the perfectly absorbing boundary condition by the method of Pedersen 11,12 , 



/ , (2^ J J i _ X (k) r . 

K M = — : : : r; 777^> C 5 - 1 ) 



d r A(fc) 



tpk 



&c(0) (27r) d J y_ w l-A(fc) 

where A(fc) is the structure factor defined by A(fc) = ^ Sf=i cos ' ^'/^) ■ For a P ar ticle starting from the first 
neighbor in Simple Cubic lattice^, 

i [ f' dd] ;?o S (k.f /b) _ 



{2ttYJ J_ v l-A(jfe) 

1 /• r m 



d d k — -±—L— = 0.516386. 



(2ir) d J J_ n 1 _ A(fc) 

Eq. (|5.1|1 together with Eq. 12.13|) for 3/> re (s) yields again the recombination probability as a function of the intrinsic 
reaction rate 7 rc at the origin. 



VI. CONCLUSIONS 



In the absence of reaction it is well known that the fractional diffusion equation is derived from the continuous 
time random walk model a 1 ' 2 ' 4 . Therefore, the continuous time random walk model can be regarded as a basis for the 
fractional diffusion equation describing sub-diffusive transport. In this paper a fractional reaction-diffusion equation 
is derived from a continuous time random walk model. The reaction term has a memory effect because the exit from 
the encounter distance which is described by the algebraic asymptotic form of the waiting time distribution of jump 
motion interferes with the reaction at the encounter distance. From the fractional reaction-diffusion equation thus 
derived the recombination probability is obtained, which depends on the intrinsic reaction rate for the sub-diffusive 
case as well as for the normal diffusion, unlike the result of Sung et al£. The theory of Sung et al. is based on the 
fractional reaction-diffusion equation where only the diffusion term has a memory effect and the reaction term does 
nol£. They obtained the recombination probability for the partially reflecting boundary condition which is equal to 
that obtained for the perfectly absorbing boundary condition 6 . An argument that is sometimes used to justify their 
result is that for the sub-diffusive case the mean residence time at a site is infinitely long, therefore any particle in a 
reactive zone undergoes reaction however small the intrinsic reaction rate may be. We have shown here that even if 
the mean residence time is infinite, each particle at a site has a finite waiting time of jump motion and a fraction of 
particles in the reaction zone will escape reaction, if the intrinsic reaction rate is finite. This is also clear from our 
simulation procedure and from the argument leading to Eq. (|3.2|l and its counterpart in the lattice model, Eq. (|5.1|l . 
To corroborate the above argument the fractional reaction-diffusion equation is derived, which has memory effects 
both in diffusion term and in reaction term. The analytical expression of the recombination probability has been 
derived from the fractional reaction-diffusion equation and coincides very well with the simulation data. The memory 
effect in the reaction term is due to the interference of the reaction at the encounter distance with the exit from the 
encounter distance which has dispersive kinetics in a sub-diffusive transport media. Although in the present paper 
the fractional reaction-diffusion equation has been derived for a specific initial condition, the linearity of the model 
guarantees a much wider range of applicability: essentially, any initial distribution will be acceptable. Finally, we 
remark that in our study the reaction is assumed to take place at some distance. Most of the theories on reactions in 
sub-diffusive transport are developed on the basis of more macroscopic description^, where the reaction is accounted 
for by the law of mass action which is space independent. In those theories, reaction terms without memory effect 
are simply added to the diffusion term with memory kernel. In this paper reactions are accounted for by the waiting 
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time distribution functions in the reaction zone and the memory effects in reaction term as well as in diffusion term 
automatically emerge from such description. 
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